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We propose a new approach to the Casimir effect based on ciassicai ray optics. We define and 
compute the contribution of classical optical paths to the Casimir force between rigid bodies. We 
reproduce the standard result for parallel plates and agree over a wide range of parameters with 
a recent numerical treatment of the sphere and plate with Dirichlet boundary conditions. Our 
approach improves upon the proximity force approximation, ft can be generalized easily to other 
geometries, other boundary conditions, to the computation of Casimir energy densities and to many 
other situations. 
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Improvements in experimental methods have rekindled 
efforts to compute Casimir forces for geometries beyond 
the classic case of parallel plates 0,0- No exact expres- 
sions are known even for simple geometries such as two 
spheres or a sphere and a plane. It is therefore interest- 
ing to consider new ways of viewing the Casimir effect 
and the approximation schemes that they motivate. In 
this Letter we present a new approach based on classi- 
cal ray optics. Our approach avoids the infinities that 
have plagued Casimir calculations. Like ray optics it is 
most accurate at short wavelengths and where diffrac- 
tion is not important. Our basic result, see eq. QBJ, is 
simple and easy to implement. It coincides with the well- 
known proximity force approximation (PFA) [3j close to 
the parallel plates limit. Recently a precise numerical 
result has been obtained for the Dirichlet Casimir en- 
ergy of a sphere of radius R separated from a plane by 
a distance ajj|. This provides us an opportunity to test 
our approximation. The results are shown in Fig. J2J). 
They give us encouragement that the optical approach 
may provide a useful tool for estimating Casimir forces 
in situations where exact calculations are not available. 
We consider a scalar field of mass m satisfying the wave 
equation, (—V 2 — k 2 )4>(x) = 0, in a domain V C M 3 
bounded by disconnected surfaces, Si, S2, ... on which it 
obeys Dirichlet (or Neumann) boundary conditions. At 
the end we comment on the generalization to conducting 
boundary conditions for the electromagnetic field. The 
Casimir energy can be written as an integral over Sp(k), 
the difference between the density of states in V and the 
density of states in vacuum. This, in turn, can be related 
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to an integral over the Green's function^ 
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where oj(k) — y/c 2 k 2 + rrfic^/K 2 and G = G — Go is the 
difference between the Green's function and the vacuum 
Green's function. We work in three dimensions although 
the generalization of our results to other dimensions is 
straightforward. We look for an approximate solution of 
the wave equation that becomes exact for infinite, planar 
surfaces. We approximate the full propagator as a sum of 
terms ascribable to different optical paths, and satisfying 
the wave equation with an error of O (l/(kR) 2 ) |g, 
where R is a typical curvature of the surface (e.g. the 
radius of the sphere in the sphere+plane situation). We 
are also neglecting diffractive contributions, arising from 
the presence of sharp boundaries, and we assume the 
absence of caustics in the integration domain. We expect 
this to give an excellent approximation for the Casimir 
energy when kR is large (here k is the dominant wave 
number in the Casimir energy integral, k ~ l/a where a 
is the minimum distance between the surfaces). 

The optical contribution to G(x',x,k) is given by the 
sum over optical paths from x to x', G(x' , x, k) — > 

^optical 

(x',x,k) = Yin G n (x', x, k). These fall into 
classes, C n which have n points on the boundaries |8|]. 
The collective index n = (n, a) identifies both the class 
C n and the path a. These paths are stationary points in 
the class C n of the functional integral representation of 
G. In three dimensions the optical terms in G contribute, 

Goptieal<y,Z,fc) 4E( 4 'V A n(^^" Ul ' l) , 
n 

(2) 

(in d 3 Hankel functions will appear in the analagous 
expression 9] ) where £ n (x',x) is the length of the opti- 
cal path n = (n, a) that starts from x and arrives at 
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x' after reflecting n times from the boundary. These 
paths are the minima of £ n (x',x), straight lines that re- 
flect with equal angles of incidence and reflection from 
the surfaces. The factor (—1)™ implements the Dirichlet 
boundary condition. For Neumann boundary conditions 
it is absent. An (a/, x) is the enlargement factor of classi- 
cal ray optics (also equal to the VanVleck determinant 
as defined in Ref. 0), and is given by 

A n (x',x) = ^ = limrt^M^+ife) , ( 3 ) 

where Ri t 2{s) are the radii of curvature of the wavefront 
following the path and s is the coordinate along the path. 
A measures the spread in area dA at the arrival point x' 
of a pencil of rays having angular width dfl at the starting 
point x, following the classical path indexed by n. This 
is reasonably easy to compute even for multiple specular 
reflections on any curved surface. The contribution of 
the optical path n to the Casimir energy is obtained by 
substituting eq. (J2J) into eq. Q, 

f°° dk f , 

£„ = (-l) n Af„ / -^hkoj(k) d 3 xJAJx) sin kUx), 
Jo 4?r Jv n 
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where M n is the multiplicity of the n-th path, and we 
have defined A n (x) = A n (x,x) and £ n (x) = £ n (x,x) for 
brevity. The minimum in the class Co (the direct path 
for x to x') should be excluded to account for the sub- 
traction of the vacuum energy. In a given geometry the 
optical paths can be indexed according to the number of 
reflections from each surface. For example in a geometry 
consisting in only two convex plates (iSi and S2) we have 
paths reflecting once on S\ or once on 1S2, paths reflect- 
ing two times (once on 1S1 and once on and so on. 
The multiplicity of even reflections paths is 2 (the path 
can be run in two different directions) while that of odd 
reflections is 1. £> n is the domain over which the path 
n = (n, a) is possible. £ n given by eq. I0J diverges if 
paths of arbitrary small length can occur. For domains 
bounded by convex plates £ n (x,x) — » can only occur 
for the first reflection, n = 1. To regulate this divergence 
we separate the initial and final points by a distance e, so 
that £ n > e. This is equivalent to putting a cutoff on the 
frequency at k ~ 1/ e fiol Hlj . Because it is confined in 
the first reflection, the divergence will never contribute 
to the force between surfaces. In practice it can be iso- 
lated and discarded. Next we interchange the integrals 
over k and x in eq. and perform the k- integral, 

£ n = / d 3 x^^-K 2 (mc£ n (x)/h). 

(5) 

For a massless scalar we let m — > and we obtain our 
fundamental result, 
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FIG. 1: a) Optical paths for parallel plates. Initial and final 
points have been separated for visibility; b) Optical paths for 
plane + sphere. For n = 1 a reflection off the sphere is shown; 
for n = 3 a reflection twice off the sphere is shown. 



which expresses the optical approximation to the Casimir 
effect as a sum over geometric quantities alone. Our 
method should not be confused with Gutzwiller's semi- 
classical approximation to the density of statesfl^. Il4|. 
nor with Balian's and Bloch's multiple reflection expan- 
sion for the Green's function The latter expresses 
the Green's function in terms of surface integrals, not 
limited to the optical paths. The former corresponds 
to performing the integration over x in eq. by sta- 
tionary phase, and selects only periodic paths. This ap- 
proximation fails badly when the radius of curvature R 
of the surface(s) is large compared both to the separa- 
tion a and the width L of the surfaces. Our method 
applies also in situations in which no periodic classical 
paths exist. Parallel plates provide a simple, pedagogical 
example which has many features — fast convergence, 
trivial isolation of divergences, dominance of the even 
reflections — that occur in all the geometries we have 
analyzed. We assume for simplicity that the two plates 
have the same area S. The relevant paths are shown 
in Fig. Q where the points x and x' , which should coin- 
cide, are separated for ease of viewing. For the even paths 
^2n( z ) = 2na, n — 1,2, . . ., independent of z (here z is 
the distance from the lower surface). For the odd paths 
^2n— i,a(z) = 2(n — l)a + 2£, where £ = z, a — z respec- 
tively if a = down, up, n — 1, 2, . . . For planar boundaries 
the enlargement factor is given by A n = \j£\- The sum 
over even reflections, 

= "5 g j Sd dS I dZ j2nW - -yS^ S (7) 

is trivial because it is independent of z. The result is 
the usual Dirichlet Casimir energy 2] . The sum over odd 
reflections gives 

£odd = 7. — o / dS > / dz 7T 

277 J i^J* (e 2 + (2z + 2na) 2 ) 2 
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The divergence as e — > can be related to that expected 
on the basis of Ref. [lOj and will be discussed in [ljj • For 
now it suffices that it is proportional to S, independent 
of a, does not contribute to the Casimir force, and can 
be ignored. The fact that the odd reflections sum to a 
divergent constant is universal for geometries with planar 
boundaries, and to a good approximation is also valid 
for curved boundaries. Notice that in this situation our 
method coincides with the method of images 0] . This 
will not occur for other examples. 1 

Note that the sum over n, eq. Q), converges rapidly: 
92% of the effect comes from the first term (the two re- 
flection path) and > 98% comes from the two and four 
reflection paths. This rapid convergence persists for all 
the geometries we an analyzed due to the rapid increase 
in the length of the paths. Also notice that for to > 
the two reflection contribution gives a uniform approx- 
imation to the a dependent part of the Casimir energy 
(accurate from 92% for ma c 1 to exponentially small 
terms for ma > 1): £\ Rn itc — £ 2 = --i^TgK 2 {2mca/h). 
The optical approach sheds light on the proximity force 
approximation, which has been used for years to estimate 
Casimir forces for geometries in which an exact calcula- 
tion is unavailable For the Dirichlet problem and two 
bodies S\ and 52 it takes the form: 

fpFA(Si) = / dS—^— (9) 

1440 J Sl [di 2 (:r)] d 

where di2 (x) is the distance from Si to £2 along the nor- 
mal to 5i at x. The PFA is ambiguous because a different 
result is obtained by interchanging surfaces S\ and £2- 
The PFA can be viewed as the sum over optical paths 
if at each point x in T> the path is chosen normal to S\ 
and the small area element that intersects this path on 
S 2 is replaced by a plane normal to the path. Then all 
paths that bounce back and forth between these two par- 
allel surfaces are summed. Clearly the PFA misses three 
important effects that are correctly included in the opti- 
cal approach a) the actual optical paths are shorter; b) 
the surfaces are curved; and c) there are optical paths 
through points that do not lie on straight lines normal to 
one surface or the other. Effects a) and c) increase and 
b) decreases for convex (increases for concave) surfaces 
the optical estimate of the absolute value of the Casimir 
energy relative to the PFA. In the cases we have studied 
the net effect is to increase the Casimir energy. In the 
subsequent discussions we compare our results with the 



The same calculation was presented as an application of 
Gutzwiller's trace formula |12| to the Casimir effect in Ref. |14| . 
In general the optical paths are clos ed but not periodic, and our 
approach will not resemble Ref. Il4l in other geometries. 
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FIG. 2: Casimir energy for a sphere and a plate as a func- 
tion of £ = a/R. Numerical data normalized to unity 
(stars with error bars) ; Optical approximation (red triangles) ; 
plate based PFA (green diamonds); sphere based PFA (blue 
squares) 



PFA approximations based on either of the two surfaces. 
First, however, we propose an "optimal" PFA motivated 
by the optical approach: at each point in V choose the 
unique shortest path from Si to £2 (of length £12). Re- 
place both surfaces locally by planes perpendicular to 
this path and sum all optical contributions. The result, 

& -< s > = -S//*« < 10 ' 

resolves the ambiguity in the PFA in favor of the shortest 
paths. Of course, the sum over the actual optical paths 
including the enlargement factor is more accurate still. 
The only non-trivial geometry we know of for which we 
can test our approximation is a sphere of radius R placed 
at a distance a from an infinite plane. This has not been 
solved analytically, but numerical results have been pub- 
lished recently |4|. Defining £ = a/R, we expect the opti- 
cal approximation to give an error of O (£ 2 )- Certainly, 
when £ > 1 diffraction dominates, the force is given by 
Casimir-Polder^ll, and our optical approximation will 
fail. Our aim is to study the accuracy of the optical 
approximation and the domain of its applicability com- 
pared, for example, to the PFA. We have calculated the 
Casimir energy for this configuration including paths up 
to four reflections. Some characteristic paths are shown 
in Fig. QJ. The results are plotted in Fig. J2J). The C\ 
and C3 contributions can be evaluated analytically. The 
divergent contribution from C\ is independent of a and 
can be put aside. The a-dependent, finite part of £\ and 
£3 are opposite in sign and their sum is always small 
(< 2%) compared to the £2- £2 and £4 can be com- 
puted quickly with Mathematica© . The optical result 
agrees with the numerical result of Ref. Q within error 
bars (~ 1%) out to £ « 0.1, where the PFA fails badly. 
Even for £ rj 1, on the border of its range of validity, 
SE/E N = 25% where E N is the numerical value given by 
Ref. [4J and SE = E N — E opt . In comparison, the "sphere 
based'' PFA gives SE/E N = 58% and the "plate based" 
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FIG. 3: Local contributions to the optical Casimir energy for 
a plane and a sphere with a/R — 0.25. The scale is linear in 
the hue from red (least) to violet (greatest). The self-energy 
density given by Ci have been subtracted. 



PFA 5E/E N = 73% 2 . The limiting case £ -> is that 
of two infinite parallel plates. Here the optical approx- 
imation must agree analytically (and numerically) with 
PFA. The agreement is clear in Fig. J3J ■ Since the optical 
approximation gives the Casimir energy as a volume inte- 
gral of a local contribution at each point x, it is possible 
to get an idea of the domains that give the dominant con- 
tributions to the Casimir force by plotting the integrand 
in eq. 10. As an example we show a contour map of the 
integrand for the plate and sphere at £ = 0.25 in Fig. |J3}. 
The local Casimir energy density (and other local observ- 
ables) is define d by a differential operator acting on the 
Greens function[ll| and is not identical to the integrand 
in eq. JHl ■ However, since the form of the optical approxi- 
mation to the Greens function is so simple, it is straight- 
forward to obtain a compact and computable approxi- 
mation to the local energy density which would replace 
Fig. ©. Our work suggests many possible extensions. 
We have studied the application to a finite (rectangular) 
plane inclined at an angle 9 to an infinite plane. This is 
the geometry of a "Casimir torsion pendulum" [l7j . The 
physically interesting case of conducting boundary con- 
ditions can be realized by constructing a matrix Green 
function for the vector potential. For parallel plates the 
contribution of odd reflections integrates to (for any a) , 
reflecting the well-known absence of a 1/e 3 divergence for 
conducting boundary conditions The even reflections 
sum to obtain the expected result: £ con d = 2£Dirichiet- 
The extension to more complicated geometries will be 
presented in Ref. |l7j . There are many interesting cases 
(both from a theoretical and an experimental point of 



2 In response to this paper, the authors of Ref. Il4l have compared 
their scmiclassical approximation with the numerical results of 
Ref. and find agreement comparable to ours ll5l . This can be 
understood as a consequence of the fact that the effective width 
of the sphere, L, scales like its radius R, and is not a general 
feature of the semiclassical approach. 



view) in which diffraction effects become important. Es- 
pecially in situations in which the objects are small com- 
pared to their separation (transition between Casimir 
and Van der Waals forces). Diffraction effects are cer- 
tainly well beyond the PFA approach but can in princi- 
le be included in the our framework by using Keller's 
recipe for constructing the diffracted rays contribu- 
tion to the propagator and then integrating over k to 
leave an x integral. The Casimir integral over modes can 
be generalized to give the partition function for a fluc- 
tuating field at finite temperature. The computation is 
straightforward and gives the thermal properties of the 
Casimir force — assuming that it remains reasonable to 
idealize the material by boundary conditions throughout 
the range of temperatures of interest 01 • 
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